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ABSTRACT 



Context. The Kepler mission has recently discovered a brown dwarf companion transiting one member of the M4V+M5V visual 
binary system LHS 6343 AB with an orbital period of 12.71 days. 

Aims. The particular interest of this transiting system lies in the synchronicity between the transits of the brown dwarf C component 
and the main modulation observed in the light curve, which is assumed to be caused by rotating starspots on the A component. We 
model the activity of this star by deriving maps of the active regions that allow us to study stellar rotation and the possible interaction 
with the brown dwarf companion. 

Methods. An average transit profile was derived, and the photometric perturbations due to spots occulted during transits are removed 
to derive more precise transit parameters. We applied a maximum entropy spot model to fit the out-of-transit optical modulation as 
observed by Kepler during an uninterrupted interval of ~ 500 days. It assumes that stellar active regions consist of cool spots and 
bright faculae whose visibility is modulated by stellar rotation. 

Results. Thanks to the extended photometric time series, we refine the determination of the transit parameters and find evidence of 
spots that are occulted by the brown dwarf during its transits. The modelling of the out-of-transit light curve of LHS 6343 A reveals 
several starspots rotating with a slightly longer period than the orbital period of the brown dwarf, i.e., 13.13 + 0.02 days. No signature 
attributable to differential rotation is observed. We find evidence of a persistent active longitude on the M dwarf preceding the sub- 
companion point by ~ 100° and lasting for at least ~ 500 days. This can be relevant for understanding how magnetic interaction 
works in low-mass binary and star-planet systems. 
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1. Introduction no evidence of chromospheric activity and presents a modula- 
J"" 1 ■ tion of very low amplitude in the Kepler optical light curve (see 
l> Kepler is a photometric space mission devoted to searching for Sectfj rj The modulation remains stable for more than 
^ planets by the method of transits and asteroseismology. The re- 500 days and has a period of 1Z71±a28 daySi which is syn . 
O cently released public data contain a large amount of informa- chronized with the orbital motion of the compan i on . The large 
m *i°n useful for characterizing stellar pulsations and modulations separation between the A and c componen ts makes tidal inter- 
t-H caused by starspots in solar-like and low-mass stars. action particularly weak in mis system. Some other known syn- 
> . M-type stars are particularly interesting targets in exoplanet chronous system s are CoRoT-3 (iDeleuil et al. 2008; Tria ud et al. 
. r-j searches since Earth-like planets in their habitable zones are ex- [200I. CoRoT-4 (lLanza et al]l2009bl). an d r Bootis ( Catala et al. 
^ pected to he within the detection capabilities of some of the cur- | W alker et al. Il2008r. IDonati et all l2008h . although in the 
U rently working instruments. However, the effects of magnetic last two cases the star is orbited by a massive planetj not by a 
^ activity in these stars may sometimes mask or even mimic the brown dwarf (here after BD). Since they are F-type dwarf stars, 
signal of a transiting exoplanet (e.g., [Bonomo & Lanza ! | 2008|; h has been sugges ted that the close-in brown dwarf or the plan- 
i Barnes et al. || 201 1|) . High-precision photometry from space of- ets may have synchronized their thin outer convecti ve envelopes 
fers the opportunity to study the signature of spots and faculae during the main _sequence lifetimes of the systems (IDonati et al.1 
during large uninterrupted intervals of time, so it represents the 2008) 
best way to improve our knowledge of activity patterns in low- 
mass stars in order to further optimize exoplanet searches. In this work we study the photospheric activity of 
Among the M-type stars observed by Kepler, a particularly LHS 6343 A by applying the spot model ling approach alrea dy 
interesting object is LHS 6343 AB, a visual binary system con- used for o ther objects, such as CoRoT-4 (Lanza et al. 2009b|) or 
sisting of an M4V and an M5V st ar at a projected distance C0R0T-6 ( lLanza et al.l 1201 ll) . This allows us to map the longi- 
of ~ 40.2 AU. iJohnson et al.l d201 ll) report the discovery of a tudinal distribution and evolution of spots in a reference frame 
brown dwarf companion (hereafter LHS 6343 C) with a mass co-rotating with the companion orbital period. As we discuss in 
of 62.7 ± 2.4Mj up transiting component A of the system every Sect. 14.21 we find an overall migration of the spot pattern, indi- 



12.7138 days. The primary star, LHS 6343 A, shows almost eating a rotation period of 13.137+0.01 1 days for LHS 6343 A, 
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and a persistently active longitude, locked at a phase difference 
of ~ 100° from the sub-companion point. Starspots are sys- 
tematically enhanced when they approach and cross that longi- 
tude during their migration. This leads us to discuss possible 
magnetic interactions between LHS 6343 A and its brown dwarf 
companion (cf. Sect. O. 



2. Observations 

About 156,000 targets are being continuously monitored by the 
Kepler space telescope i n a 10° x 10° field near the constella - 
tions of Lyra and Cygnus dKoch et al.l2010tlBorucki et alj20 10). 
Photometry is continuously acquired for time intervals up to 
~ 90 days, called "quarters" in Kepler jargon, after which a re- 
orientation of the spacecraft is required to keep the solar arrays 
pointed toward the Sun. This maneuver implies a rotation of the 
field of view on the focal plane making the image of a given 
star to fall on a different CCD. As a consequence, a jump ap- 
pears in the raw time series at the end of each quarter. Kepler 
acquires photometry with a so-called long-term cadence of ~ 30 
minutes, except for a subset of selected targets for which the ex- 
posure is reduced to about one minute. Long-cadence data have 
been made publicly available as part of the first qu arters data re- 
lease (Q 0-Q6) and their reduction is described in Jenkins et al. 
(I2010allbh . 

LHS 6343 (KIC 010002261) has been observed by Kepler 
during all the first six quarters. The time series contains a total 
of 22976 data points with ~30 minute cadence and a mean rel- 
ative precision of 7 x 10~ 5 , ranging from May 2009 to Septem- 
ber 2010. Both A and B com ponents of the system, separated 
by 0."55 dJohnson et al.ll201lt) . are contained inside the Kepler 
photometric mask used to measure the flux from the target. Its 
aperture is optimized by the Kepler pipeline for each individual 
target. Note that since the image scale is 3. "98 per pixel it is not 
possible to measure separately the fluxes coming from the A and 
B components even redefining the photometric mask. 

Cotrending basis vectors were applied to the raw data using 
PyKE pipeline reduction softwareQ to correct systematic trends, 
which are mainly related to the pointing jitter of t he satellite, de - 
tector instabilities and environment variations (Mur phvll2012l) . 
These are optimized tasks to reduce Kepler Simple Aperture 
Photometry (hereafter SAP) datcQ as they account for the po- 
sition of the specific target on the detector plane to correct for 
systematics. From two to four vectors were used for each quar- 
ter to remove the main trends from the raw data (see Fig. [!}. A 
low-order (< 4) polynomial filtering was then applied to the re- 
sulting data for each quarter because some residual trends still 
remained followed by discontinuities between quarters. These 
are due to the change of the target position on the focal plane 
following each re-orientation of the spacecraft at the end of each 
quater. Raw SAP data and detrended light curves are shown in 
Fig-Q] As a consequence of this data reduction process, the gen- 
eral trends disappeared, as well as jumps of the photometry be- 
tween quarters. On the other hand, any possible intrinsic long- 
term variability of the object (> 50 days) has also been removed. 
The use of low-order polynomials to detrend each quarter series 
of data ensures that the frequency and amplitude of any vari- 
ability with a timescale comparable with the companion orbital 
period is preserved. The presence of a number of gaps in the data 




1 http://keplergo.arc.nasa.gov/PyKE.shtml 

2 The SAP light curve is a pixel summation time-series of all calibrated 
flux falling within the optimal aperture. 
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Fig. 2. Lomb-Scargle periodogram of the out-of-transit light modula- 
tion of LHS 6343. A Gaussian fit was applied to the highest peak to ob- 
tain the period of the main modulation and its uncertainty (12.71 ± 0.28 
days). Peaks corresponding to the harmonics of the main periodicity 
also appear at higher frequences. 



prevents us from using other detrending methods such as Fourier 
filtering. 

The out-of-transit light curve was cleaned by applying a 
moving box-car median 3cr-clipping filter to identify and remove 
outliers. Then it was rebinned with a regular 120 min sampling 
to reduce the computation load for the spot modelling. This is al- 
lowed because variations related to spot activity have timescales 
of the order of one day or longer. The average standard error 
of the binned observations is ~ 5.0 x 10~ 5 in rela tive flux units. 
Transi ts were removed assuming the ephemeris of [Johnson et alJ 
d201 ll) . The resulting light curve consists of 5716 data points 
ranging from BJD 2454953.569 to BJD 2455462.163, i.e., with 
a duration of 508.593 days. 

Studying the photometric variability of LHS 6343, we note 
that the light curve exhibits a periodic modulation with a low 
amplitude (~ 10~ 3 in relative flux units, resulting in a signal-to- 
noise ratio around 15 - 20) that remains stable throughout all 
the time series. The minima of the modulation are separated 
by a nearly constant phase interval from the transits of the BD 
companion, as can be seen in the lower panel of Fig. Q] which 
also shows the long-term coherence of the modulation. This 
suggests that component A is responsible for the observed vari- 
ability. A first analysis of the out-of-transit light curve was per- 
formed through a Lomb-Scargle periodogram, finding only one 
significant period of 12.71 + 0.28 days, thus confirming the syn- 
chronicity with the LHS 6343 C orbital period (see Fig. |2). The 
other peaks seen in the periodogram correspond to the harmon- 
ics of the main periodicity and have heights similar to that of the 
main peak because the signal is quasi-periodic, but not exactly 
sinusoidal in shape. A sinusoidal fit to the light curve with the 
above period was performed to estimate the phase lag, showing 
that the light minima precede the mid-transits by 98° in phase. 
The phased light curve obtained with all the Q0 - Q6 data is pre- 
sented in Fig [3] with overplotted a sinusoidal best fit, showing 
the difference in phase between the minima of the photometric 
modulation and the transits. 

For further analysis and spot modelling, we assume that all 
the observed variability is caused by LHS 6343 A. This is also 
supported by the fact that the system is estimated to be older 
than 1-2 Gyr from the very low observed chromospheric activity, 
and most likely near 5 Gyr from brown dwarf mass and radius 
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Fig. 1. Upper panel: Kepler SAP light curve of LHS 6343, covering the first six quarters (red dots), and after correction from main trends using 
cotrending basis vectors (grey dots). Lower panel: LHS 6343 light curve after low-order polynomial filtering, as adopted for subsequent analysis. 
The times of mid-transits of LHS 6343 C are marked by vertical dashed lines. 



relations (iJohnson et al.ll201 U iBaraffe et alj|1998l) . Therefore, 
the rotation period of the isolated component LHS 6343 B is ex- 
pected to be longer than that of the observed modulation (iBarnesI 
l2003ll2007h . 

Since both A and B components of the system fall within 
the Kepler photometric aperture, it is necessary to correct for the 
flux contribution from star B before pe rforming any modelling 
on the light curve. IJohnson et al.l (1201 ll) use the measured B — V 
colours and V magnitudes to estimate the relative magnitudes of 
the stars in the Kepler bandpass, obtaining AKp = 0.74 ±0.10. 
This yields a flux ratio of 1.97 ± 0.19 between star A and B. 
By assuming that all the modulation comes from star A, we use 
this value to correct for the dilution of the flux produced by the 
component B. This hypothesis is valid if the amplitude AFb of 
the variability of component B satisfies the relationship: 



AF R 



AF A 



Fa + Fn 



<3C 



Fa + Fn 



(1) 



where AFa is the variability of component A and Fa + Fb the 
total flux of the system. Given the above flux ratio, this implies: 



AF B AF A 
— ^« (1.97 ±0.19)— ± 



(2) 



In other words, the relative amplitude of the variability of 
LHS 6343 B must be much lower than twice that of the A com- 
ponent to ensure that our analysis is consistent. This seems to be 
reasonable in view of the motivations given above. 



3. Light curve modelling 



3.1. Transit modelling 



Although lJohnson et al.l d201 ll) present a transit light curve anal- 
ysis, recomputating the system parameters considering the cur- 
rently available dataset is justified, since the more extended Q0- 
Q6 transit photometry reveals several features related to spot oc- 
cultations that would cause a bias in the computation of the tran- 
sit depth if not properly accounted for. Specifically, when the 
transiting BD occults a dark feature on the disc of the A compo- 
nent, a relative flux increase appears along the ligh t curve le ading 
to a re duced relative depth of the transit (cf., e.g., iBallerini et al.l 
I2012I) . 

For modelling the transits, we retained the original Ke- 
pler time cadence of ~ 29.4 minutes without applying the 3cr- 
clipping filter to the time series. Our fi tting procedure for the 
phased transit photometry is based on the Mandel ^fc Agoll (f2002) 
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Fig. 3. Upper panel: Kepler light curve of LHS 6343 phased with 
a period of 12.71 days (black dots). The transits have been removed 
leaving a gap centred on phase 0. The sinusoidal best fit applied to 
calculate the phase lead (~ 98°) is plotted as a red solid line. Lower 
panel: The vertical scale is zoomed in and the Kepler data are plotted 
in 0.025 phase bins to highlight the deviations from a sinusoidal shape 
which cause the appearance of harmonics of the main frequency in the 
periodogram of Fig. [2] Error bars are 1-cr deviations of the binned flux 
values. 



analytic mod el and is performe d by using the Transit Analy- 
sis Package dGazak et al J 120 1 2h . a graphic user-interface IDL 
tool. This combines a Markov chain Monte Carlo (MCMC) tech- 
nique to fit light curve s and a wavelet-based likelihood function 
dCarter&Winnll2009h . which is used to compute parameter un- 
certainties. The model light curve is resampled and rebinned 
according to the Kepler long-cadence (~ 29.4 minutes) data 
sampling to avoid systema tic deviations, i n particular during the 
ingress and egress phases (Kipping 2010). The orbital period is 
fixed to 12.71382 days (IJohnson et alJl201lh . and the adjusted 
parameters are the inclination i, the ratio of the companion and 
primary radii Rq/Ra, the scaled semimajor axis (o/Ra), and the 
time of mid-transit 7q. Limb darkening is treated by adopting 
a quadratic approximation, and the coefficients are interpolated 
for the specifi c LHS 6343 A stellar pa rameters from the tabu- 
lated values of IClaret & Bloerrienl (1201 ll) . An attempt to derive 
the limb-darkening coefficients from the data fitting while fix- 
ing the rest of the parameters has resulted in unphysical values, 
probably due to a complex brightness profile and the presence of 
active regions. Both u\ and ui have therefore been fixed to the 
values obtained from the tables. In most cases, the low impact of 
the limb darkening on the shape of the transit - note its almost 
flat-bottomed profile - and the existing degeneracy with other 
parameters, such as the inclination, makes this a better approach 
than obtaining the coefficients from transit modelling. 

The eccentricity and argument of periastron are given by 
IJohnson et al.l d201 ll) . A low value of the eccentricity (e 0.05) 
may still affect the transit modelling by increasing its duration by 
~ 5% when the argu ment of periastro n is changed from to — 0° 
to to = 180° (cf. lPerrvmanll20TTl) . Indeed, there is a cer- 
tain probability that the eccentricity is zero because the value 
is within two standard deviations from zero (assuming a Gaus- 
sian distribution of the parameter). Moreover, the eccentricity is 
a definite positive quantity so any orbital fit to a radial velocity 
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Fig. 4. Upper panel: Kepler LHS 6343 transit photometry from Q0- 
Q6 (black dots) and the best fit model obtained as described in Sect. 13.11 
(red line); crosses indicate data corresponding to spot occultations dur- 
ing the transits. They are not used for modelling the transit. Lower 
panel: The residuals of the transit best fit (black dots) and the points 
corresponding to spot occultations (crosses) vs. the time elapsed from 
mid-transit. 

Table 1. Parameters obtained from the transit light curve modelling of 
LHS 6343 A. 



Parameter 


Value" 


Radius ratio (Rc/Ra) 


0.221 ± 0.002 


Semimajor axis (a,R = 0,/Ra) 


45.5 + 0.4 


Orbit inclination (i) 


89.60 + 0.04 


Impact parameter (b) 


0.30 ± 0.03 


Mid transit (BJD-2454957.0) 


0.2174 + 0.0001 


Linear limb darkening (u^a) 


0.457* 


Quadratic limb darkening (m2,a) 


0.374* 


Eccentricity 


0.056 ± 0.032 c 


Argument of periastron, to (deg) 


-23 ± 56 f 


Mean stellar density, pa (p©) 


6.8 + 0.4 



Uncertainties correspond to 85.1% confidence interval from 
MCMC fitting plus the contrib ution of systematic error s. 
b Fixed values obtained from lClaret & Bloemenl(l201 ll) . 
c Fixed values obtained from IJohnson et all (1201 ll) . 

time series with some degree of noise tends to give a non-zero 
value for this parameter ev en if the orbit is perfectly circular (cf., 
e.g., iHusnoo et al. Il2012h . The eccentricity is expected to be 
zero in our system, which is close to synchronization, but a pos- 
sible excitation by perturbations caused by the distant B com- 
panion cannot be excluded in the absence of a specific study. 
However, a transit analysis with the assumption of a circular or- 
bit was also performed and obtained no significant change in the 
parameter values. 

A first modelling of the transit was performed by using a 
Levenberg-Marquardt algorithm, and a three-sigma clipping fil- 
ter was applied to the positive residuals of this fit to identify 
and remove the data points affected by occulted spots. Then, an 
MCMC fitting algorithm was applied and the resulting posterior 
parameter probability distribution was obtained from the combi- 
nation of ten chains. Several extra runs were performed by fixing 
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each of the parameters in turn to check the consistency of the re- 
sults and find possible systematic errors that were added to the 
final uncertainties. The error on the stellar flux ratio (see Sect.[2]i 
was also accounted for in the presented uncertainties. The pres- 
ence of surface brightness inhomogeneities, associated with stel- 
lar activity, makes it especially difficult to accurately model the 
transit with our long-cadence data because the induced pertur- 
bations are often unresolved. This produces an increase in the 
in-transit residuals. The results are presented in Table [T]and the 
fit displayed in Fig. [4] The posterior probability density func- 
tions of the parameters obtained by the MCMC modelling are 
shown in Fig.0assuming uniform priors. 

Our mean transit profile, derived from 41 observed transits of 
LHS 6343 C after removing the spot perturbations, has a better 
sampling of the ingress and egress intervals allowing us to re- 
duce the uncertainties on th e parameters in comp arison with the 
previous determination by Johnson et al.1 Q2 11). Specifically, 
we find a sligthly lower value for the relative radius and a higher 
value for the inclination, although both are within two standard 
deviations of their previous values, respectively. That we find 
a very similar result for Rc/Ra essentially means that spot oc- 
cultations did not occur or at least were not sign ificant for the 
five tra nsit events present in the data analysed by Johnso net al.1 
d201 ll) . On the other hand, the results for the relative orbital 
semimajor axis a/R\ are in close agreement. Also the mean 
stellar density as computed fr om the transit parameters (e.g., 
Seager & Malle n-Ornelasll2003l) is consistent with the mass and 
radius parameters in Table|2] adopted from Johns on et al.l(l20 1 ll) . 

Some systematic trends are observed in the residuals, al- 
most symmetrically placed around the second and third contacts. 
They could be due to a complex brightness profile along the oc- 
culted chord on the stellar disc. In other words, a quadratic limb- 
darkening law could not be a good approximation for this very 
late-type star. Short-cadence data should be used to investigate 
these systematic deviations further along individual transits. 

The occultation of the stellar disc by the BD companion 
during the transits can be used to resolve the fine structure of 
the spot pattern and study its ev olution (e.g., ISilva-Valio et al.l 
2010: ISilva-Valio & Lanzall201 ll) . From the phase-folded data 
shown in Fig. [4] we find evidence of occulted spots that seem 
to primarily appear close to mid-transit. This could mean 
that some active regions are associated with the sub-companion 
point. Therefore, short-cadence data for LHS 6343 would be 
particularly interesting for further study. This would also help 
for constraining the projected obliquit y and confirm the quasi- 
synchronicity of the system (c f., e.g., ISanchfs-Oi eda et ap 201 1: 
ISanchfs-Oieda & WinrJ20Tlh . 

Transit residuals are analysed to look for variability in the 
transit depth. One of the effects producing such a variability 
could be a change in the third light contributed by the compo- 
nent B. We find that transit residuals are dominated by noise and 
do not show any modulation or trend, thus confirming that star 
B is not active enough to produce a detectable variability and is 
probably rotating with a much longer period than star A. How- 
ever, since some scatter in the transit depth may be caused also 
by active regions in star A that are not occulted during transit, 
it is difficult to use this result to put a strong constraint on the 
variability of LHS 6343 B. 

3.2. Spot modelling 

The reconstruction of the surface brightness distribution from 
the rotational modulation of the stellar flux is an ill-posed prob- 
lem, because the variation in the flux vs. rotational phase con- 



tains information only on the distribution of the brightness in- 
homogeneities vs. longitude. The integration over the stellar 
disc effectively cancels any latitudinal information, particularly 
when the inclination of the rotation axis along the line of sight 
is clo se to 90°, as is ass umed in the present case (see Sect. 13.31 
and lLanza etaf 2009a). Therefore, we need to include a pri- 
ori information in the light curve inversion process to obtain a 
unique and stable map. This is done by computing a maximum 
entropy (hereafter ME) map, which has been proven to success- 
fully reproduce active re gion distributions and area variations in 
the case of the Sun (cf. lLanza et al. 20071) . 

Our appro ach has been comp ared with other spo t mode lling 
methods by iMosser et all (120091) and iHuber et all (120101) and 
with the maps obtained by spot occ ultations during the planetar y 
transits in the CoRoT-2 system by ISilva-Valio & Lanzal d201 lb . 
To model the high-precision space-bome photometry of e Eri 
and k l Ceti obtained by MOST (the M i crovariability and O scil- 
lation of STar satel lite), ICroll I d2006b . ICroU et al.l (120061) . and 
IWalker et al.l d2007l) have developed a sophisticated approach 
based on a few circular spots. They use simulated annealing and 
tempering techniques to explore the multi-dimensional model 
parameter space and Markov chain Monte Carlo methods to as- 
sess the parameter uncertainties. The main limitation of these 
models is that the reduced x 1 values of their best fits are sig- 
nificantly greater than the unity. This indicates that the model 
assumptions, in particular the consideration of only two or three 
fixed spots, are too simplified to reproduce the true stellar vari- 
ability down to the level of precision of space-bome photometry. 

A similar a pproach has been ap plie d in the case of the Kep ler 
light curves bv lFrasca et"aTI (1201 ll) and lFrohlich et al.l (1201 2|) by 
considering several spots whose areas evolve with time. Al- 
though these models allow a convincing determination of the 
model parameters and their uncertainties, they require a huge 
amount of computation and their best fits often show system- 
atic deviations suggesting that the spot pattern is more com- 
plex than can be described with several circular, non-overlapping 
starspots. Therefore, we prefer to use a continuos distribution for 
the spotted area and a regularization approach to find a unique 
and stable solution. 

An important advantage of the ME regularization is that the 
multi-dimensional parameter space of the objective function to 
be minimized (cf. Eq. [6} has a smoother landscape than the x 1 
landscape of a multi-spot model, thanks to the addition of the 
configurational entropy term (the function S defined below; cf. 
the an alogous problem of image restoration in iBrvan &"Sk illing 
1980). This effectively reduces the problems associated with the 
presence of several x 1 relative minima falli ng in widely sep- 
arated regions of the parameter space (cf. ICroll et all [2006 ; 
IWalker et al.l2007l) . Moreover, we effectively average over many 
degrees of freedom by only considering the distribution of the 
spotted area vs. the stellar longitude, which is the only informa- 
tion that can be directly extracted from our unidimensional time 
series. Therefore, the basic role of our regularization is that of 
smoothing the distribution of the spot filling factor versus lon- 
gitude, thus reducing the impact of the noise on the light curve 
inversion process. 

In our model, the star is subdivided into 200 surface ele- 
ments, namely 200 squares of side 18°, with each element con- 
taining unperturbed photosphere, dark spots, and facular areas. 
The fraction of an element covered by dark spots is indicated by 
the filling factor /, the fractional area of the faculae is Qf, and 
the fractional area of the unperturbed photosphere is 1 - (Q + 1 )/. 
The contribution to the stellar flux coming from the A:-th surface 
element at the time tj, where j = 1, JV is an index numbering 
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Fig. 5. Probability density functions for the transit parameters obtained with the MCMC algorithm. The dotted vertical lines correspond to 2 - cr 
confidence levels. 



the N points along the light curve, is given by 

AF kj = Io(Mkj){l-(Q + l)fk + c s fk+ 

QfkW + c f (l - fikj)]}AklXkjwQi kj ), (3) 

where Iq is the specific intensity in the continuum of the unper- 
turbed photosphere at the isophotal wavelength of the observa- 
tions; Cs_jnd_£f_are_the spot and facular contrasts, respectively 
(cf. lLanza et"a"fl l2()04): is the area of the k-th surface ele- 
ment, 

, . ( 1 if fi kj > ... 
w ^ = \0 if^<0 (4) 

is its visibility; and 

\i k j = cos iffy = sin i sin 6k cos[4 + Q.(tj - to)] + cos i cos 9k, (5) 

is the cosine of the angle iffy between the normal to the sur- 
face element and the direction of the observer, with i being the 
inclination of the stellar rotation axis along the line of sight; 
9k the colatitude and 4 the longitude of the A:-th surface el- 
ement; Q. denotes the angular velocity of rotation of the star 
(Q = 2n/P ml , with P rot the stellar rotation period), and to the ini- 
tial time. The specific intensity in the continuum v aries accord- 
ing to a quadratic limb-darkening law, as adopted by lLanza et alj 
(l2003h for the case of the Sun, viz. Iq cc a p + b p fj. + c p fi 2 . The 
stellar flux computed at the time tj is then F(tf) = £t AF^. To 
warrant a relative precision of the order of 10~ 5 in the compu- 
tation of the flux F, each surface element is further subdivided 
into 1° x l°-elements, and their contributions, calculated accord- 
ing to Eq. (0, are summed up at each given time to compute 



the contribution of the 18° x 18° surface element to which they 
belong. 

We fit the light curve by varying the value of / over the sur- 
face of the star, while Q is held constant. Even fixing the rotation 
period, the inclin ation, and the spot and facular contrasts (see 
lLanza et"a"fll2007l for details), the model has 200 free param- 
eters and suffers from non-uniqueness and instability. To find a 
unique and stable spot map, we apply ME regularization by min- 
imizing a functional Z, which is a linear combination of the x 1 
and the entropy functional S ; i.e., 

Z=^ 2 (f)-^(f), (6) 

where f is the vector of the filling factors of the surface elements, 
A > a Lagrangian multiplier determining the trade-off between 
light cu rve fitting and re gularization. The expression for S is 
given in lLanza et al.l dl998h . The entropy functional S attains its 
maximum value when the star is free of spots. Therefore, by in- 
creasing the Lagrangian multiplier A, we increase the weight of 
S in the model, and the area of the spots is progressively reduced. 
This gives rise to systematically negative residuals between the 
observations and the best-fit model when A > 0. The optimal 
value of A depends on the information content of the light curve, 
which in turn depends on the ratio of the amplitude of its rota- 
tional modulation to the average error of its points. In the case 
of LHS 6343 A, the amplitude of the rotational modulation is 
~ 0.001, while the average error of the points is ~ 5 x 10 in 
relative flux units, giving a signal-to-noise ratio of ~ 20. The 
significant impact of the noise on the light curve makes the reg- 
ularization a critical proce ss, and a diff erent criterion than that 
applied for CoRo T-2 (see lLanza et aL 2009a) or CoRoT-6 (see 
lLanza etaT]l201 lh must be adopted, because in those cases the 
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signal-to-noise ratio was > 100. For the present case, we in- 
crease the regularization until we obtain |// reg | = /feo, where ^ reg 
is the mean of the residuals of the regularized light curve, and 
e = cr u / y/N is the standard error of the residuals of the un- 
regularized light curve (i.e., obtained with A = 0), defined as 
the ratio of their standard deviation <r to the square root of the 
number of data points N in each individual segment of duration 
Aff fitted by our model (see below). The parameter (3 > 1 de- 
pends on the information content of the light curve and is unity 
when the signal-to-noise ratio is of the order of 100 or greater. 
In our case, we determine the optimal value of f3 a posteriori, 
by increasing A until the small random spots associated with the 
noise are strongly smoothed, while the deviation of the regular- 
ized model from the data points is still acceptable (see Sect. 14.11 
and Appendix [A] for the effects of varying /3 on the solution). 

In the case of the Sun, by assuming a fixed distribution of the 
filling factor, it is possible to obtain a good fit of the irradiance 
changes only for a limited time interval Aff, not exceeding 14 
days, which is the lifetime of the largest sunspot groups domi- 
nating the irradiance variation. In the case of other active stars, 
the value of Aff must be determined from the observations them- 
selves, looking for the longest time interval that allows a good 
fit with the applied model (see Sect. 13.31 ). 

The optimal values of the spot and facular contrasts and of 
the facular-to-spotted area ratio Q in stellar active regions are 
unknown a priori. In our model the facular contrast Cf and the 
parameter Q enter as the product Cf Q, so we can fix Cf and vary 
Q, estimating its best value by minimizing the x 2 of the model, 
as shown in Sect. 13.31 Since there are many free parameters in 
the ME model is large, for this spec ific application we make use 
of the model of Lanza et al. (2003), which fits the light curve 
by assuming only three active regions to model the rotational 
modulation of the flux plus a uniformly distributed background 
to account for the variations in the mean light level. This pro- 
cedure is the same as used to fix the value of Q in the cases of 
CoRoT -2, CoRoT-4, CoRoT-7, and CoRoT-6 (cf. lLanza et alj 
l2009allbl 2010. 2011). 

We assume an inclination of the rotation axis of LHS 6343 A 
of i = 90° in most of our models (see Sect. l3.31 >. coming from the 
results of the transit modelling (see Sect. I3.lt . Such a high in- 
clination implies that negligible information on the spot latitude 
can be extracted from the rotational modulation of the flux, so 
that the ME regularization virtually puts all the spots around the 
subobserver latitude (i.e., w 0°) to minimize their area and max- 
imize the entropy. Therefore, we are limited to mapping only 
the distribution of the active regions vs. longi tude, which can be 
done w ith a resolution better than ~ 50° (cf. lLanza etaT 2007, 
2009a). Our ignorance of the true facular contribution may lead 
to systematic errors in the active region longitudes derived by 
our model because faculae produce an increase in the flux when 
they are close to the limb, leading to a systematic shift of the 
longitudes of the active regions used to reproduce t he observed 
flux modulation, as discussed by Lanza et alj ( 2007b for the case 
of the Sun and illustrated by Lanz a et al. (l2009al cf. Figs. 4 and 
5) for CoRoT-2. 

Given the dependence of our modelling approach on several 
parameters, we present a detailed discussion of their impact on 
our results in Appendix lAl 

3.3. Model parameters 



Table 2. Parameters adopted for the light curve modelling of 
LHS 6343 A. 



Parameter 




Ker. 


otar mass \JV1q) 


U. J /U 


T1 1 
J 1 1 


Cf v *.n Jinn / D \ 

Mar radius (/<0j 


U.J 10 


J 1 1 


I eff l^V) 




T 1 1 
J 1 1 


logg (cm s~ 2 ) 


4.85 


J 1 1 




r\ i /:o.£T 

u.ioyj 


n i i 
LBl 1 


bp 


1 OA/1 1 

1.2041 


rr) 1 1 
Col 1 


Cp 


-0 3736 


V I_> 1 1 


P mt (days) 


12.7138 


Jll 


e 


6.09 x 10" 6 


Hel2 


Inclination (deg) 


89.60 


Hel2 


Cs 


0.536 


Hel2 


Cf 


0.115 


L04 


Q 


0.0, 8.0 


Hel2 


Aff (days) 


8.769 


Hel2 



a Referenc es: CB1 1: IClaret & Bloemenl (1201 ll); Hel2: p r 
work; Jll: iJohnson et alJ (1201 ll) : L04: lLanza et all <|2004). 



resent 



The fu ndamental stellar parameters are taken from lJohnson et all 
(1201 ll) and are listed in Table|2] A quadratic limb-darkening law 



is adopted for the stellar photosphere, viz., oc « p + £> p /j + 
Cp// 2 , where the coefficients a p , b p , and c p have been obtained for 
T e fi = 3130K, \ogg = 4.85 [cgs ], and solar abundances fro m the 
theoretical values computed by IClaret & Bloemenl (1201 lb from 
Phoenix model atmospheres. 

The rotation period for modelling the spotted photosphere 
was fixed to the companion orbital period. This assumption 
comes from the periodogram analysis giving P rot = 12.71 ± 
0.28 days (see Fig. [2), which indicates that the system is close to 
synchronization. As we see in Sect. 14.21 any deviation of stellar 
rotation from synchronization will appear on the ME map as a 
longitude drift of the active regions versus time. 

The rotation period of 12.71 days is used to compute the po- 
lar flattening of the star due to the centrifugal potential in the 
Roche approximation, yielding a relative difference between the 
equatorial and the polar radii e of 6.09 x 10~ 6 . The effect of this 
flattening on the modulation caused by spots is negligible. The 
relative difference in flux for a ~ 2% spot coverage would be 
less than 10~ 6 as a consequence of the gravity darkening in the 
equatorial regions of the star. 

Owing to the lack of a measurement of the Rossiter- 
McLaughlin effect or precise v sin i data, the inclination of the 
stellar axis of LHS 6343 A cannot be constrained . The limit on 
v sin i < 2 km s _I given by IJohnson et alj d201 lb does not help 
because the above rotation period and the radius yield an equato- 
rial rotation velocity of ~ 1.5 km s _1 . Nevertheless, we assume 
that the system is spin-orbit aligned; that is, its inclination is as- 
sumed to be equal to that of the orbit of the BD companion, i.e. 
89?60 + ?04. The inclination of the stellar spin axis along the 
line of sight is the only parameter required to characterize the 
geometry of the stellar rotation for modelling the out-of-transit 
light modulation. If we had been interested in also modelling 
the occultations of starspots during transits, the sky-projected 
angle between the stellar spin and the orbi tal angular momen- 
tum w ould have been required, too (e.g., ISanch is-Oieda et al. 
120121) . 

The maximum time interval Aff that our model can accu- 
rately fit with a fixed distribution of active regions has been de- 
termined by means of a series of tests dividing the total interval, 
T = 508.5936 days, into N f equal segments, i.e., Af f = T/Nf. 
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Fig. 6. The ratio of the x~ of the composite best fit of the entire 
time series to its minimum value vs. the parameter Q, i.e., the ratio of 
the facular area to the cool spot area in active regions. The horizontal 
dashed line indicates the 95% confidence level for y 2 /y 2 . This does 
not put any significant constraint on the interval of acceptable Q values 
owing to the low signal-to-noise ratio of the data. 

The simple three-spot model of lLanza et al.l (120031) has been ap- 
plied to fit the individual intervals to reduce the number of free 
parameters. The optimal Nf is given by the minimum value that 
allows us a good fit to the light curve as measured by the x 1 
statistics. For Nf < 58, the quality of the fit starts to degrade 
significantly, so we adopt this as the optimal number of seg- 
ments, giving Aff = 8.769 days (cf. Appendix lAl. In any case, 
Atf > frot/2 is required to adequately sample the light modu- 
lation produced by a rotating active region, providing us with 
information on its longitude and area. On the other hand, a too 
long Aff would result in a low time resolution of our model, so 
that starspot evolution and migration would be more difficult to 
trace. 

The spot tem perature contrast is estimated by the correlation 
found in Fig. 7ofBerdyug]ni(|2005j) that shows that the contrast 
decreases with decreasing stellar effective temperature. Thus, we 
adopt a spot temperature deficit of AT S = -250 K for LHS 6343 
A and compute the spot's intensity contrast considering black- 
body emission for the spots and the photosphere in the Kepler 
passband, yielding c s = 0.536. This is not a critical parameter 
because a variation in c s will only cause a change in the over- 
all spot filling factor, but will not affect either spot longitudes 
or their time evolution, as discussed in detail by lLanza et all 
(2009a|). T he facular contras t is assumed to be solar-like with 
c f = 0.115 (lLanza et all2004l> . 

The ratio Q between the facular and the spotted areas in ac- 
tive regions is estimated by m eans of the three-spot model in- 
troduced bv lLanza et al.l (120031) simultaneously with the optimal 
value of Aff. Both the flux increase due to the facular compo- 
nent when an active region is close to the limb and the flux de- 
crease due to its spots when that region transits across the central 
meridian of the disc are modelled by this approach, allowing us 
to derive the optimal value of Q, if the signal-to-noise ratio of 
the data is high enough. 

In Fig. [6] we plot the ratio X 2 /Xmm °^ tne total ^ 2 of the com- 
posite best fit of the entire time series to its minimum value x^ n , 
versus Q. Given the low signal-to-noise ratio of the photometry 
in our case, we find that there is no clear ^-minimum vs. Q so 
its value is poorly constrained with the best solutions falling in 
the range between ~ 5 and ~ 11. Therefore, we adopt Q = 8.0 
for further analysis. This facular-to-spotted area ratio is similar 



to the one found in the case of the Sun (Q = 9, cf. Lanz a et aD 
l200l . 

4. Results 

4.1. Light curve model 

The composite best fit obtained with the ME regularization and 
Q = 8.0 is shown in Fig. |7](upper panel) together with the resid- 
uals (lower panel). From a best fit without regularization (A = 0) 
applied to the whole time series (Q0-Q6), we obtain a mean 
// res = 8.308 x 10~ 8 and a standard deviation of the residuals 
(To = 6.778 x 10~ 5 in relative flux units, which is slightly higher 
than the standard error of the mean points. This can be due to 
short time scale fluctuations that cannot be fitted with our ap- 
proach that reproduces variations on a timescale comparable to 
Aff or longer. On average, we have N ~ 98 data points per fitted 
subset Aff. 

The Lagrangian multipliers A for the regularized ME mod- 
els are iteratively adjusted until the mean of the residuals jj. reg = 

-1.365 x 10" 5 =i -y3(o- / VP), where ft = 2 is established as 
the optimal value that still leads to an acceptable fit while re- 
ducing the information content of the spot maps as much as 
possible, i.e., maximizing the entropy. A lower value of f3 im- 
proves the fit slightly over small sections of the light curve, but 
introduces several small spots that change from one time inter- 
val to the next without any clear regularity. On the other hand, 
by selecting /? = 2, we obtain smoother maps with concentrated 
groups of spots that evolve regularly from one Aff interval to the 
next, which indicates that we are properly modelling the over- 
all distribution of the active regions on the star. A higher level 
of regularization (J3 > 2) gives unacceptably larger, systemati- 
cally negative residuals because the maximum entropy criterion 
reduces the spotted area too much. At the same time, it con- 
centrates the spots into isolated groups, and a pattern with three 
active regions equally spaced in longitude (i.e, with w 120° gaps) 
begins to appear. This pattern still shows the general migration 
of the multiple spots appearing in the maps with the optimal reg- 
ularization (see Sect. 14.2b . but it is clearly an artefact of the over- 
regularization, due to the limited signal-to-noise level of our data 
that favours a reproduction of the light variations by a decreasing 
number of spots when the regularization is increased (cf. Ap- 
pendixlAl. 

In either case, the low signal-to-noise ratio of our data re- 
quires a relatively strong regularization, and as a consequence, 
misfits sometimes appear at the matching points between time 
intervals. Therefore, by decreasing Aff, the overall quality of 
the best fit is not increased in our case. Nevertheless, look- 
ing at the residuals of our regularized composite best fit, we 
see a generally good reproduction of the light variations with- 
out systematics that could be attributed to sizable active regions 
evolving on timescales significantly shorter than Aff. We notice 
that the noise level increases after BJD 2455276 (correspond- 
ing to BJD-2454833 = 443 in Fig. |7j, which coincides with 
the beginning of quarter Q5 and returns to the previous level af- 
ter BJD 2455371 (BJD-2454833 = 538), i.e., at the beginning 
of quarter Q6. This is a consequence of the rotation of the fo- 
cal plane of the Kepler telescope between successive quarters, 
which makes the target be measured on a different CCD. 

4.2. Distribution of active regions and stellar rotation 

The distribution of the spot filling factor / versus longitude and 
time is plotted in Fig. [8]for our ME model with Q — 8.0 adopting 
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Fig. 7. Upper panel: The out-of-transit light curve of LHS 6343 A (dots) and its ME-regularized best fit for a facular-to-spotted area ratio 
Q = 8.0 (solid line) during the time interval between BJD 2454953.569 and BJD 2455462.163. The flux is normalized to the maximum observed 
flux along the time series. Lower panel: The corresponding residuals. 



a reference frame rotating with the orbital period of the BD com- 
panion, i.e., 12.7138 days. The origin of the longitude frame cor- 
responds to the subobserver point at the beginning of the Kepler 
time series (i.e. BJD 2454953.569) and the longitude increases 
in the same direction as the stellar rotation and the orbital mo- 
tion of the companion. The subcompanion longitude is fixed at 
~ 257° in this reference frame. 

Several individual active regions can be identified from the 
map in Fig. [8]showing lifetimes around 15-20 days. Three dot- 
ted lines have been drawn connecting the longitudes with higher 
filling factors to trace the overall migration of the spots. They 
identify three long-lived active longitudes that migrate back- 
wards at the same constant rate. This indicates that the rotation 
of the star is somewhat slower than that of our reference frame 
following the orbital motion of the BD companion. The simi- 
larity of the migration rates among the three active longitudes 
suggests that either the spot latitudes are confined to a thin belt 
or the star has no detectable differential rotation. 

By considering the longitudes of the relative maxima of the 
spot covering factor in Fig. [8] a linear fit is performed to find 
the relative rotation rates associated with the steady migration 
of the three main starspots. The results are listed in Table [3] 
for two models having a different facular-to-spotted area ratio 
Q. By assuming a significant facular contribution (Q — 8.0), a 
mean rotation period of 13.137 ± 0.01 1 days is obtained. This is 
a significantly slower rotation than the companion orbital mo- 
tion, showing that the system is slightly out of synchronism. 
The stellar rotation, as traced by starspots, appears to be rigid 
since there is no significant difference between the highest and 
the lowest migration rates. The migration rates for a spot model 
assuming only dark spots (Q - 0) are also reported in Table [3] 
because the longitudes of the spots derived from th e maximum 
entropy mod elling may depend on the value of Q (lLanza et alj 
2007, 2009a). Nevertheless, very similar results are obtained, 
thus confirming that the presence of faculae is not actually re- 



Table 3. Relative migration rates AQ./Q obtained from linear fits to the 
migration of the main active longitudes in Fig. [8] and the corresponding 
mean rotation period P mt . 



Q 


AL] 


AL 2 


AL 3 






(10~ 2 deg/d) 


(10~ 2 deg/d) 


(10~ 2 deg/d) 


(d) 


0.0 


-3.20 ±0.20 


-3.34 ±0.11 


-3.24 ±0.17 


13.129 ±0.015 


8.0 


-3.34 ±0.15 


-3.34 ±0.14 


-3.30 ± 0.20 


13.137 ±0.011 



quired to adequately fit the low signal-to-noise photometry of 
LHS 6343 A (cf. Sect. l3.3b . In the case of no faculae, a virtually 
identical mean rotation period of 13. 129 ±0.015 days is obtained 
from the migration of the spots. 

An intriguing phenomenon seen in Fig.[8]is that the spot fill- 
ing factor shows a persistent relative minimum in the longitude 
range between * 100° and m 250°, while it has a relative max- 
imum between w 250° and w 450°. In other words, the filling 
factor appears to be systematically enhanced when the migrating 
starspots approach and cross a longitude leading the subcom- 
panion longitude by ~ 100° in our reference frame. Since the 
passage of that longitude across the centre of the stellar disc pre- 
cedes that of the BD companion, this agrees with the light curve 
that shows a relative light minimum ~ 3.5 days before each tran- 
sit. This modulation has persisted for more than 500 days, as far 
as the Kepler data currently cover. 

This enhancement of the spot filling factor is highlighted in 
Fig- El where the spot coverage factor integrated over the time is 
plotted vs. longitude in bins of 72°. A wide maximum is evident 
around ~ 350°, meaning that activity tends to be concentrated 
in a region preceding the subcompanion longitude. This agrees 
closely with what can be seen in Fig. Q~| and especially in the 
phase-folded light curve (Fig. |3). Moreover, this explains the 
results obtained by the Lomb-Scargle periodogram (Fig. |2) be- 
cause only the frequency of the main modulation associated with 
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Fig. 8. Isocontours of the ratio f/f mllx , where / is the spot covering 
factor and / m;lx = 0.00059 its maximum value, versus time and longi- 
tude for the ME models with Q = 8.0. The two dashed black lines mark 
longitudes 0° and 360° beyond which the distributions are repeated to 
easily follow spot migration. The dashed white line marks the longi- 
tude of the subcompanion point, which is fixed in this reference frame 
at ~ 257°. The contour levels are separated by A/ = 0. l/ milx with yel- 
low and orange indicating the maximum covering factor and dark blue 
the minimum. The dotted white lines trace the migration of the active 
regions associated with each active longitude (see the text). 



this activity enhancement, phased to the companion C, is de- 
tected by that method. On the other hand, ME modelling allows 
us to map individual starspots and follow their migration across 
the longitude where they show an enhancement. Unfortunately, 
no information on the spot latitudes can be derived with our mod- 
els, so we cannot establish whether the enhancement of activity 
at this longitude is produced in a single latitude or if it covers a 
range of latitudes. As discussed in Sect. [5] these results strongly 
suggest some magnetic interaction between LHS 6343 A and its 
BD companion LHS 6343 C, thus causing spots to mainly appear 
at a fixed longitude in the reference frame of the companion or- 
bital motion. 

The variation in the total spotted area vs. the time is shown 
in Fig. [10] Only slight variations can be seen with a typical 
timescale of 20-30 days, which can be associated with the typi- 
cal spot lifetime. Short-term oscillations are most likely related 
to the noise present in our time series, especially during quarter 
Q5, residual trends, or gaps in the data (e.g., the one near day 
~ 400, see Fig. [TJ. Long-term trends have been removed when 
detrending the time series from Kepler instrumental effects (see 
Sect. 13, so the mean value of the total spotted area stays ap- 
proximately constant. The absolute value of the area depends on 
the adopted spot contrast c s and the value of Q (for instance, a 
lower spot temperature would imply a stronger contrast and thus 
a smaller area), but the relati ve variations in the area are largely 
independent of c s and Q (cf. lLanza et"a l. 2009a). 
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Fig. 9. The spotted area in 72° longitude bins averaged versus time 
as derived from the ME regularized models with Q = 8.0. The error 
bars indicate the standard error of the spotted area in the corresponding 
bin. The dotted line marks the subcompanion point longitude, while 
the dashed lines indicate the longitudes 0° and 360° beyond which the 
distribution is repeated for clarity. 
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Fig. 10. The total spotted area versus time as derived from the reg- 
ularized ME analysis for Q = 8.0. The error bars correspond to three 
standard deviations. 

5. Discussion 

We find that the rotation of LHS 6343 A is not perfectly syn- 
chronized with the orbital motion of the BD companion. Spot 
modelling is superior to the Lomb-Scargle periodogram for de- 
termining the rotation period because a much smaller phase dif- 
ference is required to distinguish between close periods for a 
given time baseline. The regular and steady migration of the 
active longitudes in Fig.[8]makes our result quite robust. Unfor- 
tunately, the fast evolution of the individual spots together with 
the low signal-to-noise ratio of the photometry does not allow us 
to detect any possible differential rotation. 

The intriguing point concerning the quasi-synchronous rota- 
tion of LHS 6343 A is that the tidal interaction between the star 
and its BD companion is pre sently very weak. Using the formal- 
ism of lLeconte et al. I (12010b . we estimate a current timescale for 
the tidal angular momentum exchange between the spin and the 
orbit of as 7.5 Gyr, if a modified tidal quality factor Q' = 10 6 is 
adopted for the A component. Assuming that the stellar radius 
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was larger by a factor of 3 - 4 during the pre-main-sequence evo- 
lution of the system, the synchronization timescale was shorter 
by a factor of » 10 3 at that time, i.e., significantly shorter than 
the pre-main-sequence lifetime of the A component. Therefore, 
the system should have reached the ZAMS close to synchroniza- 
tion and with the stellar spin aligned with the orbital angular 
momentum. Then it evolved to the present state with a likely 
modest angular momentum loss from the A component. Indeed, 
a reduced angular momentum loss rate in comparison with ear- 
lier main-sequenc e stars seems to be a characteristic of mid-M 
dwarfs (cf., e.g., iReiners et al.l 12012b . Alternatively, the pres- 
ence of a close-in companion may have reduced the efficiency 
of the magne tic braking by the stellar wind, as conjectured by 
Lanza (2010) in the case of stars accompanied by hot Jupiters 
(cf. also ICohen etal. 112010b . 

The possibility that magnetic interactions in star-planet sys- 
tems lead to an enhancement of st ellar activity was initially 
suggested bv lCuntz etal] (l2000h and lCuntz & Shkolnikl §002). 
Some evidence has been found in several star-planet systems 
from the modulation of chromospheric activity indexes in phase 
with the or bit of the planet rath er than with the rotation period 
of the star (IShkolnik et al.ll2005l) . For HD 179949, HD 189733, 
and t Boo, an enhancement of the flux variability has been ob- 
served preceding the s ubplanetary point by ~ 70° during several 
observing campaigns (IShkolnik et al.ll2008l) . 

In the case of LHS 6343 A, where the companion is a low- 
mass BD object, the available data are photometric, so the ac- 
tivity modulation is observed in the photosphere. Here, a phase 
lead of 98° between the transits and the photometric minima is 
observed (see Sects. [2] and I4.2l i. The effect peaks once per or- 
bit, thus suggesting a magnetic interaction because two relative 
maxima per orbit are expected in the case of a tidal interaction 
due to the symmetry of the tidal bulge. Moreover, the relative 
amplitude of the flux modulation due to the ellipsoidal deforma- 
tion of LHS 6343 A does not exceed ~ 5 x 10~ 5 owing to th e 
large separation of the A and C components dPfahl et al. | |2008). 
Even considering an orbit with an eccentricity of 0.1, the differ- 
ence in the ellipsoidal light variation between the periastron and 
the apoastron does not exceed ~ 3 x 10~ 5 in relative flux units. 
Therefore, the observed flux modulation with a relative ampli- 
tude of ~ 10~ 3 , phased to the orbital motion of the BD, cannot 
be attributed to the tidal distortion of the primary star. 

Similar cases of photospheric active regions rotating in phase 
with a close-in mas sive planet hav e been reported fo r t Boo 
dWalker et al. I f2008) and CoRoT-4 (lLanza et al.ll2009bb . In the 
former case, a photospheric active region producing a signature 
of the order of one millimag was observed in different seasons 
with the MOST satellite. The feature was sometimes dark, some- 
times bright and was associated in phase with the persistent chro- 
mospheric hot spot discussed above. In the case of CoRoT-4, 
a preferential longitude for spot appearance was found within 
~ 40° - 50° from the subcompanion point. The rotation of both 
stars is close to synchronization as in the case of LHS 6343 A. 
Although t Boo and CoRoT-4 are F-type stars, the same mecha- 
nism could be at work in all three cases. 

Several models have been proposed to account for the fea- 
tures of the magnetic interaction between stars and close-in plan- 
ets. They can be extended to the case o f a close BD compan - 
ion as in LHS 6343 AC. The models b y iMcIvor et al. I (12006). 
iPreusse et al.l (120061) . and lLanzal d2008l) have focussed on inter- 
preting of the phase lag between the planet and the maximum 
visibility of the chromospheric hot spot by adopting different 
hypotheses, e.g., a tilted potential stellar magnetic field, a fi- 
nite propagation time of the magnetic perturbation along the field 



lines from the planet to the star, or a non-potential twisted field 
in the stellar corona. They are generally capable of account- 
ing for the phase lag observed in the case of a chromospheric 
hot spot, but cannot account for an enhancement of the photo- 
sph eric activity. As sh own by the numerical MHD simulations 
of (ICohen et all 120091 2011a,b), when the planetary magneto- 
sphere is inside or in touch with the stellar Alfven surface, a re- 
markable interaction between the stellar and the planetary fields 
occurs, and the energy released in the magnetic reconnection 
can be transferred down to the stellar chromosphere and pho- 
tosphere. This can account for a localized enhancement of the 
flux, i.e., a local brightening as in large solar flares that have a 
white-light photospheric counterpart to their chromospheric and 
coronal emission. However, the appearance of dark spots associ- 
ate d with a close-by com panion requires a different mechanism. 

lLanzal (120081 1201 ll) have conjectured that the perturba- 
tion of the coronal field and the triggering of flaring activity 
in the corona may increase the loss rate of magnetic helicity 
through coronal mass ejections (cf. also lLanz al l2012l) leading 
to a perturbation of the stellar dynamo. As reviewed by, e.g., 
iBrandenburg et all (1201 2|) . the operation of the dynamo in late- 
type stars implies a mechanism that can get rid of the magnetic 
helicity that is continuously pumped into the small field scales, 
otherwise the dynamo is rapidly quenched by magnetohydrody- 
namic effects. In the Sun, coronal mass ejections associated with 
large flares are regarded as a viable mechanism for accomplish- 
ing this goal, therefore an enhancement of flaring activity in stars 
with close-in companions may help the operation of the dynamo. 
Since the perturbation moves in phase with the planet, the associ- 
ated loss of magnetic helicity is non-axisymmetric and may lead 
to a modulation of the a effect of the stellar dynamo in p hase 
with the motion of the companion (see lLanzal 12008, 201 li for 
details). In systems that are close to synchronization, the mod- 
ulation can be particularly effective thanks to the slow relative 
motion of the companion in the reference frame rotating with the 
star that increases the duration of the interaction between a stel- 
lar coronal structure and the magnetic field of the orbiting body. 
The phase lag between the starspots and the companion may be 
explained by adopting a twisted magnetic fie ld interconne cting 
the spots with the orbiting BD as discussed in Lanza (2008J). 

An alternative explanation taking into account the destabi- 
lizing effects of the tidal perturbation on the magnetic flux tubes 
inside the red dwarf star cannot of course be excluded. In this 
scenario, the persistent active longitude would correspond to the 
periastron passage of the BD companion with the phase lag of 
~ 98° being the angular distance of that point from the mid- 
transit. Recent works suggesting a possible gravitational per- 
turbation of the planets on th e solar dynamo m ay give support 
to this conjecture (cf., e.g., lAbreu et aO2012[ and references 
therein ). 

6. Conclusions 

Considering 508 days of Kepler photometry of the visual system 
LHS 6343, we have computed a new best fit to the mean pho- 
tometric transit profile of the C brown dwarf component across 
the disc of the A comp onent. We essentially confirm the tran- 
sit parameters found by [Johnson et al.l (1201 ll) and reduce their 
uncertainties. 

A maximum entropy spot model was applied to the Kepler 
light curve of LHS 6343 to derive maps of the longitudinal distri- 
bution of the active regions on the A component, assumed to be 
responsible for the out-of-transit modulation of the flux. While 
including faculae is not found to be critical for obtaining an ac- 
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ceptable fit to the light modulation, an optimal fit of the light 
curve is found for a facular-to-spotted area ratio around Q — 8.0. 

Our maps show several active regions with a slower rota- 
tion rate (P rot = 13.14 ± 0.02 days) than the BD orbital pe- 
riod (P or b = 12.71 days), implying that the AC system is quasi- 
synchronized. Owing to the limited signal-to-noise ratio of the 
data and the fast evolution of the active regions that have life- 
times of 1520 days, it is not possible to detect any differential 
rotation signature. On the other hand, the starspot covering fac- 
tor is observed to be systematically enhanced when the spots 
cross a longitude leading the substellar point by ~ 100°. This 
effect is also detected in the photometric time series as a modu- 
lation with the orbital period of the BD showing minima ~ 3.5 
days before each transit. That this active longitude persists for 
more than ~ 500 days provides evidence that the BD companion 
does affect the photospheric activity of its host somewhat. We 
conjecture that a magnetic interaction between the coronal field 
of the M dwarf and the BD may produce a modulation of the 
stellar dynamo leading to an increase in activity modulated in 
phase with the orbital motion (see Sect. O. 

A monitoring of the system through measurements of the 
modulation of the chromospheric lines (e.g., Ha or Ca II H&K) 
and coronal X-ray emission would be very helpful, because a 
detection of a hot spot in the upper stellar atmosphere with the 
same phase lag as the detected photospheric active longitude 
would indicate interactions through magnetic field reconnection, 
thus supporting the conjectured mechanism. 
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Appendix A: Influence of the adopted model 
parameters on the spot maps 

Appendix A.1: Time resolution (At f ) and amount of 
regularization (J3) 

In Sect. I3.3I we discussed the method applied to determining 
the maximum time interval Aff that our model can accurately 
fit when considerin g a fixed distribution of three active regions 
(lLanza et alj|20 03). This corresponds to the time resolution of 
our model, since we assume that active regions do not evolve on 
timescales shorter than Aff . The total interval T was divided to 
Nf segments so that Aff = T/Nf . The optimal Nt was found to be 
58 from the best fit of the three-spot model, as measured from 
the^f 2 statistics. 

Three cases corresponding to different lengths of the time 
interval Aff are presented in Fig. IA.1I where the ratio X 2 /Xmm 
is plotted vs. Q, the facular-to-spotted area ratio; x 2 is me c hi 
square of the composite best fit to the total time series; and^- 2 nin 
is the minimum^ 2 obtained for Nj = 58 and Q = 8.0. The cor- 
responding spot model is discussed in the text. Here we explore 
the case with Nf = 50 and perform an ME analysis with the cor- 
responding time interval, i.e., Aff = 10.1719 days, to investigate 
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Fig. A.l. The ratio of the^ 2 of the three-spot model fit to the minimum 
value for Nf = 58 vs. the parameter Q, i.e., the ratio of the facular 
area to the cool spot area in active regions. Solution series for three 
different values of the time interval Aff = T/Nf are shown. 

the effect of a different Aff on our results. The same approach as 
described in Sect. 14. II was considered, and regularized ME mod- 
els were iteratively adjusted until /3 = 2. The isocontour map 
of the distribution of the spot filling factor vs. time and longi- 
tude is presented in Fig. IA.2l It is remarkably similar to the one 
presented in Fig. [8] and the phased activity enhancement around 
~ 350° is still visible, although the migration of the starspots and 
their changes on short timescales are less evident. Therefore, a 
longer time interval limits the possibility of accurately measur- 
ing spots' migration and estimate their lifetimes. 

In a second test, we chose the same time interval as in Sect. [4] 
i.e., Nf = 58, to explore the effects of a different amount of 
regularization on the ME spot maps. In Fig. IA.3I we present a 
spot map obtained with a lower regularization (/3 — 1). At that 
point, the large starspots appearing in Fig. [8] are resolved into 
several smaller spots. This fine structure changes vs. the time 
without any clear regularity, as expected for an artefact due to 
overfitting the noise present in the data. 

Finally, we plot the distribution of the filling factor obtained 
with /3 = 3 in Fig. IA.4l The effects of the overegularization de- 
scribed in Sect. 14.1 l ean be seen as a tendency for the appearance 
of larger and smoother spot groups that include previously re- 
solved starspots. Three main groups separated by approximately 
120° in longitude generally appear, showing the same backwards 
migration as with (3 — 2, the optimal value adopted for the reg- 
ularization. The residuals of the composite best fit to the light 
curve (not shown) show larger systematic deviations than in the 
cases with /? = 1 or/3 — 2, again a sign of overegularization. 

Appendix A.2: Inclination of the stellar spin axis 

The inclination of the stellar spin is assumed to be equal to that 
of the orbit of the BD companion in our analysis, i.e. i = 89?60. 
An isotro pic orientation of th e stellar spin axis has its mode at 
i = 90° (iHerreroet al.ll2012l) . so it is reasonable to choose a 
value close to this, given that there is no a priori information to 
constrain the parameter. 

The tidal timescale for the alignment of the stellar spin and 
the orbital angular momentum are comparable, in the case of 
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Fig. A.2. As in Fig. [8] but adopting a longer time interval Aff = 
10.1719 days. 

LHS 6343 A, to the synchronization timescale that is estimated 
to be « 7.5 Gyr. Given a probable age of the system between 1 
and 5 Gyr, we cannot be sure that the obliquity has been damped 
by tides, unless we assume that most of the damping occurred 
during the pre-main-sequence phase when the tidal interaction 
was stronger due to a larger (cf. Sect. [5}. Therefore, we 
adopt a prudential approach and explore different cases in which 
the system is out of spin-orbit alignment. 

We adopt the same parameters as in Sect. 14.1 1 and compute 
a regularized ME best fit of the light curve with /3 = 2, fixing 
the stellar inclination to i - 60° and i = 45°. The resulting spot 
maps are plotted in Figs. IA.5l and lA.6l respectively. We may ex- 
pect that the latitudinal distribution of the spots are remarkably 
different, especially for i = 45°. On the other hand, the plot- 
ted distributions of the filling factor vs. longitude and time are 
similar to what is computed with i = 89?60. Also the evolution 
and migration of the active regions are similar to that case. Only 
some details appear to be critically dependent on the adopted 
inclination. 

The results of the tests described in this Appendix confirms 
that by averaging the ME spot maps over latitude, we can effec- 
tively remove most of the degeneracy present in the light curve 
inversion process. 
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Fig. A.3. As in Fig. [8] but performing an ME regularization with 
/>' I- 
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Fig. A.4. As in Fig. [8] but performing an ME regularization with 
P = 3. 



Fig. A.5. As in Fig. [8] but adopting a stellar inclination of ( = 60°. 
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Fig. A.6. As in Fig. [8] but adopting a stellar inclination of i = 45°. 
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